Non- Fermi-liquid behavior in cubic phase BaRuC>3: A dynamical 

mean-field study 



Li Huang 1 ' 2 and Bingyun Ao 2 

1 Beijing National Laboratory for Condensed Matter Physics, 
and Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China 
2 Science and Technology on Surface Physics and Chemistry Laboratory, 
P.O. Box 718-35, Mianyang 621907, Sichuan, China 
(Dated: March 8, 2013) 

Abstract 

Motivated by the recently synthesized cubic phase BaRu03 under high pressure and high tem- 
perature, a thorough study has been conducted on its temperature-dependent electronic properties 
by using the state-of-the-art ab inito computing framework of density functional theory combined 
with dynamical mean-field theory. At ambient condition the cubic phase BaRu03 should be a 
weakly correlated Hund's metal with local magnetic moment. The spin-spin correlation function 
and local magnetic susceptibility can be well described by the Curie- Weiss law over a wide temper- 
ature range. The calculated low-frequency self-energy functions of Ru-4(i states apparently deviate 
from the behaviors predicted by Landau Fermi-liquid theory. Beyond that, the low- frequency 
optical conductivity can be fitted to a power-law 5Rcr(w) ~ u; -0 ' 98 , which further confirms the 
Non- Fermi-liquid metallic state. 

PACS numbers: 71.30.+h, 71.27.+a, 75.30.Wx 
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I. INTRODUCTION 



Recently the alkaline-earth ruthenium oxides have attracted growing interest. These 
oxides generally exhibit fascinating physics properties, such as unconventional superconduc- 
tivity with p-wave symmetry (S^RuO^,- antiferromagnetic Mott insulator (Ca2Ru04), 2 
and orbital selective Mott transition in S^-zCa^RuO^- etc. Among the rest, due to their 
interesting magnetic phase diagrams, transport properties and potential device applications, 
the ternary ruthenates with perovskite or perovskite-related structures (ARUO3: A = Ca, 
Sr, Ba) have been extensively studied by numerous experiments and theoretical calculations 
in the past decade.- - — 

Both CaRuC>3 and SrRuC>3 crystallize in the orthorhombic perovskite structure with a 
GdFeC>3-type distortion. SrRuC>3 is a highly correlated, narrow-band metallic ferromagnet 
with a Curie temperature (T c ) of about 160 K.- Its local magnetic moment (1.4 fis) is rather 
large, despite highly extended 4d character of the valence electrons. Interestingly, CaRu03, 
an isostructural compound, does not show any magnetic ordering in finite temperatures.- 
The nature of its magnetic ground state still remains controversial. We note that one of 
the most striking properties of CaRu03 and SrRu03 compounds is the violation of Landau 
Fermi-liquid (LFL) theory, which has been proven by many experimental results, including x- 
ray photoemission spectra, transport and optical properties, etc.-"- The strength of Coulomb 
interaction U and the importance of Hund's rule coupling J among Ru-4<i orbitals are 
another two interesting topics and in lively debate. Though almost all the experimental^"— 
and theoretical"— efforts manifest some role of electron-electron correlation, the strength 
and the extent of its importance still remain unclear. 

In the earlier years, it is well known that depending on how BaRuOs is synthesized it 
has several polytype structures, i.e., the nine-layered rhombohedral (9R), the four layered 
hexagonal (4H), and the six layered hexagonal (6H).— ^ Lately, the cubic phase BaRu03 
with ideal perovskite structure has been synthesized by Jin et al. under 18 GPa at 1000 
°C— It remains metallic down to 4 K and occurs a ferromagnetic transition at T c = 60 K,— 
which is significantly lower than that of SrRuOs.- The ferromagnetic transition in SrRu03 
falls into the mean-field universality class whereas cubic phase BaRuOs exhibits significant 
critical fluctuations as described by the 3D Heisenberg model.— The availability of cubic 
phase BaRu0 3 not only completes the polymorph of BaRu0 3 , but also makes it possible to 

2 



map out the evolution of magnetism and other properties as a function of the ionic size of 
the A-size in the whole series of ARUO3.— 

Despite tremendous efforts have been made, little is known about the basic properties of 
cubic phase BaRuCV In this paper, we will address the following two issues: (i) Definitely, 
in CaRu03 and SrRuC>3, the effects induced by electronic correlation can not be ignored.—"— 
However, does the electronic correlation play a decisive role in cubic phase BaRuC>3? (ii) It 
is believed that the physical properties of CaRu0 3 and SrRu0 3 can not be well described by 
the LFL theory. Indeed, the evidences are their low-frequency optic conductivity, resistivity, 
and electronic Raman scattering intensity which obey the fractional power-law.-"- Thus, 
whether the physical properties of cubic phase BaRuC>3 still violate the LFL theory is an 
essential problem. 

The density functional theory within local density approximation combined with dy- 
namical mean-field theory (dubbed LDA + DMFT) is a very powerful computing frame- 
work for strongly correlated materials.—"— In the present works, by employing the LDA 
+ DMFT computational scheme, the temperature-dependent electronic properties of cubic 
phase BaRuC>3 have been systematically studied. In contrast to CaRuC>3 and SrRuC>3, un- 
der room temperature cubic phase BaRuC>3 is a weakly correlated Hund's metal with local 
magnetic moment, and its low-frequency conductivity deviates the u~ 2 law as is predicted 
by classic LFL theory. 

II. METHOD 

We first compute the ground state electronic structures of cubic phase BaRuC>3 within 
nonmagnetic configuration by utilizing the plane-wave pseudopotential approach as is im- 
plemented in the quantum espresso software package.— The generalized gradient approxi- 
mation with Perdew-Burke-Ernzerhof exchange-correlation functional is used to describe 
the exchange and correlation potentials. The pseudopotentials in projector augmented wave 
scheme^ for Ba, Ru, and O species are built by ourselves. The electronic wave functions are 
described with a plane-wave basis truncated at 80 Ha, and a T-centered 12 x 12 x 12 fc-point 
grid for Brillouin zone integrations is adopted. These pseudopotentials and computational 
parameters are carefully checked and tuned to ensure the numerical convergences. 

To include the effect of electronic correlation, the ground state wave functions are used 
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to construct a basis of maximally localized Wannier functions (MLWF) for Ru-4g? and 0-2p 
orbitals. The corresponding multiband Hubbard Hamiltonian has the following form2i^ 



where n im = c im Ci m , and c im (ci m ) creates (destroys) an electron in a Wannier orbital m at 
site i. Here i^LDA is the effective low-energy Hamiltonian in the basis of Ru-4<i and 0-2p 
MLWF orbitals, and thus is a 14 x 14 matrix. i^DC is a double counting correction term 
which accounts for the electronic correlation already described by the LDA part, and the 
fully local limit scheme^ is chosen. The Coulomb interaction is taken into considerations 
merely among the Ru-4<i orbitals. We use U = 4.0 eV and J = 0.65 eV to parameterize 
the Coulomb interaction matrix, which are close to previous estimations.— 1 ^ To solve the 
many-body Hamiltonian (JTJ) , in the DMFT part^ 1 ^ 1 . we employ the hybridization expansion 
continuous time quantum Monte Carlo impurity solver (abbreviated CT-HYB).— »2£ Finally, 
through the mature analytical continuation methods^ 1 ^ the impurity spectral function can 
be extracted directly from the imaginary-time Green's function which are derived from the 
quantum Monte Carlo simulations. 

III. RESULTS AND DISCUSSION 

Figure [1] represents our calculated orbital-resolved density of states for Ru-4<i states at 
several typical temperatures. The octahedral surrounding of Ru splits the Ru-4g? states into 
three-fold degenerated t 2g and two-fold degenerated e g levels. Ru 4+ ion, which is nominally 
in a low-spin, d 4 configuration, gives rise to a t\ g configuration with Fermi level lying in the 
t 2g manifold with empty e g states. As for the density of states of t 2g states, it displays a 
sharp quasiparticle peak near the Fermi level, a shoulder structure around -0.3 eV, and a 
Hubbard subband like hump at -8.0 eV ~ -2.0 eV. While for the density of states of e g states, 
since it is less occupied, the primary spectral weight is above the Fermi level. There are two 
small satellites located on both sides of the Fermi level (-0.1 and 0.3 eV, respectively). With 
the increment of temperature, the two peaks will be smeared out gradually. To sum up, the 
integrated spectral functions of Ru-4<i states show significantly metallic features, and the 
temperature effect is not very obvious. When the temperature rises from 190 K to 580 K, 
slightly spectral weight transfer to high energy is observed. 
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FIG. 1. (Color online) Spectral functions of cubic perovskite BaRu03 calculated by LDA + DMFT 
method. Upper panel: Spectral functions of Ru-4<i states at various temperatures. Lower panel: 
Quasiparticle band structure of BaRu03 along high symmetry lines in the Brillouin zone. 



In the next step, we computed the full momentum-resolved spectral function A(k, u) 
along some high symmetry lines in the Brillouin zone for cubic phase BaRuOa. The inverse 
temperature (3 is chosen to be 40, which corresponding to T = 290 K approximately. In 
the lower panel of FigfTJ A(k, u) is shown in comparison with the LDA band structure. 
The sharp quasiparticle peak observed in the integrated spectral function is clearly visible 
on the intensity plot, and fairly well defined. It lies in the region about from -2.0 eV to 
1.5 eV, dominated by the t 2g states. At higher energy, the e g states become the majority. 
However, in the region below -2.0 eV, the 0-2p states make a major contribution. From 
the distribution of spectral weights of t 2g and e g states, it is speculated that in the region 
from -2.0 eV to -7.0 eV, there exists strong hybridization between the Ru-4o? and 0-2p 
states. Comparing this with the LDA band structure, first of all we notice the quasiparticle 
band structure does not show apparent shifting. Secondly, except for becoming diffuse, 
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FIG. 2. (Color online) Magnetic properties of cubic perovskite BaRu03. (a) Spin-spin correlation 
functions x( T ) = (S z (0)S z (t)) at various temperatures, (b) Inverse local magnetic susceptibility 
Xi oc as a function of temperature. 



the band renormalization of the quasiparticle band structure is hard to be distinguished. In 
general, the quasiparticle band structure of cubic phase BaRu03 coincides with its LDA band 
structure, giving rise to a picture of weakly correlated metal. On the contrary, previous LDA 
+ DMFT calculations for CaRuOa and SrRuOa present strongly renormalized and shifted 
quasiparticle band structure,— resulting in the picture of moderately correlated metal. 

In recent years, the evolutional trend of ferromagnetism in ARuC"3 is in hot debate. 
Thus in the present works, we calculated the spin-spin correlation function x( r ) an d local 
magnetic susceptibility xioc of cubic phase BaRu0 3 , and tried to elucidate its magnetic prop- 
erties in finite temperatures. The calculated spin-spin correlation functions are illustrated 
in FigfS^a). On one hand, the cubic phase BaRu03 exhibits a well-defined frozen local mo- 
ment, which is characterized by a spin-spin correlation function that approaches non-zero 
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FIG. 3. (Color online) Electronic self-energy function of Ru-4d states, (a) Imaginary part of the 
Matsubara self-energy function SjS(iw) for t2 g and e g orbitals at T = 290 K. The solid lines denote 
as the fitted function — S£(iu;) = A(iu) a + B. (b) Orbital-resolved low-energy scattering rate 
7m = -9S m (iw -» 0). 

constants at large enough t, as is easily seen from T = 193 K to 1160 K. On the other hand, 
the spin-spin correlation function does not behave as x( r ) ~ {T/ sin(Tr7r)) 2 for times r suf- 
ficiently far from r = or (3 respectively, which means the violation of LFL theory.— From 
the spin-spin correlation function, the local magnetic susceptibility Xioc = Jq x( r )dT can be 
easily evaluated, which is plotted in FigJ2^b). As shown, the calculated xioc is Curie- Weiss 
like over a rather wide temperature range, in other words, it follows a xToci^) — T/C law 
at least up to T = 1160 K. This implies that the Ru-4d electrons in cubic phase BaRuOs 
retain the local nature of the magnetic moment. 

Next we concentrate our attentions to the electronic self-energy functions of Ru-4<i states. 
The calculated orbital-resolved ^£(20;) are shown in Figj3](a). For the sake of simplicity, 
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FIG. 4. (Color online) Real part of optical conductivity of cubic perovskite BaRuOs by LDA + 
DMFT calculations. Inset: The low- frequency ^Ra(uj) at T = 290 K, and the dashed line represents 
the fitted function ?R.a(oj) = Cuj~ a . 

only those results calculated at T = 290 K are presented. Werner et al— have suggested 
that the still-mysterious optical conductivity cr(u) in pseudocubic SrRuOs and CaRuOa, 
which varies approximately as u~ ' 5 and deviates sharply from the prediction of LFL theory, 
can be perfectly interpreted as a consequence of square-root self-energy function. Inspired 
by their works, we conducted a careful analysis to determine the asymptotically formula for 
the low-frequency self-energy function. In a Fermi-liquid, the imaginary part of Matsubara 
self-energy should exhibit a linear regime at low energy, whose slope is directly related to 
the quasiparticle mass enhancement. However, as shown in Figj3]^a), we do not observe any 
linear behavior: the Matsubara self-energy behaves as — QY>(iu) = A(iu) a + 7 with a ~ 0.48 
for t 2g states and a ~ 0.80 for e g states, respectively. The non-linear frequency dependence 
of the Matsubara self-energy implies that Landau quasiparticles and effective masses can 
not be properly defined for cubic phase BaRuOs. The non-zero intercept 7 = —QIl(iuj — > 0) 
can be viewed as the low-energy scattering rate and it is a broadly used physical quantity to 
distinguish the LFL and NFL phases.— As a byproduct, the orbital-resolved 7 m is evaluated 
as a function of temperature and shown in Figj3]^b). Clearly, 7 m increases monotonously 
with the increment of temperature and ^ t2g > 7e s is always valid. For both t 2g and e g states 
7 m can not be neglected even at T = 190 K. Thus, it means that similar to SrRuOs and 
CaRu0 3 , the cubic phase BaRu0 3 lies in the NFL regime as well. 

A power-law analysis on the transport properties, such as resistivity p{T) cx T n , of cubic 
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phase BaRuOa was made by Zhou et al. 19 and the exponent n as a function of pressure 
was evaluated recently. Their results show an interesting evolution from n ~ 1.85, which 
is close to n = 2 for the LFL phase at ambient pressure, to n ~ 1.4 of the NFL phase 
at the pressure where the ferromagnetic phase collapses. The most important evidence 
for NFL state in SrRuC>3 and CaRuOs is the fractional power-law conduct ivity.- - - In this 
work, we also calculate the optical conductivity cr(u>) of cubic phase BaRuC>3 under various 
temperatures. In FigJUonly the real part of optical conductivity is shown. The sharp peak 
near u> = denotes the Drude-like feature. The broad hump located from 1.5 eV to 3.5 eV 
can be attributed to the contribution of electron transition between quasiparticle peak and 
Hubbard subbands.— With the increment of temperature, this hump slightly shifts to higher 
frequency region, which is in accord with the variation trend of Hubbard subbands observed 
in the temperature-dependent integrated spectral functions of Ru-4g? states (see FigJT]). In 
order to further confirm whether the underlying physics of cubic phase BaRuC>3 can be 
described with LFL theory, we conduct a detailed power-law analysis for the low-frequency 
optical conductivity at T = 290 K. The low-frequency optical conductivity is fitted by the 
exponent function §ta(u) ~ Cu~ a . The quantitative results are shown in the inset of FigJH 
The fitted exponent a ~ 0.98, while the expected value predicted by LFL theory is a = 2. 
It's worth mentioning that the exponent a is approximately 0.5 for pseudocubic SrRuC>3 and 
CaRu0 3 , and 0.7 for some high temperature superconductivity cuprates.-^ Nevertheless, the 
optical conductivity data suggest the NFL metallic nature of cubic phase BaRuC>3 under 
ambient condition again. 

Finally, we should emphasize the importance of Hund's physics in cubic phase BaRuC^. 
Very recent investigations about iron pnictides and chalcogenides showed that strong cor- 
relation is not always caused by the Hubbard interaction U, but can arise from the Hund's 
rule coupling J.— ^ Since the strength of electronic correlation in these materials is al- 
most entirely due to the Hund's rule coupling, they are dubbed Hund's metals by Haule 
et alM at first. It has recently been noticed by Yin et al. that in realistic Hund's met- 
als, the electronic self-energy and corresponding optical conductivity show NFL power-law 
frequency dependence, tendency towards strong orbital differentiation, and that large mass 
enhancement can occur even though no clear Hubbard subband exist in the single particle 
spectra.— 1 ^"— According to their investigations, both the iron pnictides and chalcogenides 
are typical Hund's metals. The origin of fractional power-law in the optical conductivity of 
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them can be traced to the Hund's rule interaction. As for the cubic phase BaRu03, based 
on our calculated results: NFL behavior in low-frequency self-energy function and scattering 
rate, fractional power-law in the optical conductivity, and considerable mass enhancement 
(at T = 290 K m* 2g = 1.8mo an d fn* eg = 1.2mo), we can conclude that it is another realistic 
Hund's metal. Indeed, we have performed additional LDA + DMFT calculations for cubic 
phase BaRuOs with different Coulomb interaction strengths from U = 2.0 eV to 6.0 eV 
and obtained almost identical results. However, when U = 4.0 eV and the Hund's rule 
coupling term is completely ignored (J = 0.0 eV), the NFL behaviors previously found in 
the self-energy function and optical conductivity are absent totally. It should be noted that 
Medici et a/.— have suggested that the physical properties of ruthenates are governed by the 
Hund's physics, in other words, the pervoskite ARUO3 forms a new series of Hund's metal. 
Our calculated results for the cubic phase BaRuC>3 confirm their issue as well. 

IV. SUMMARY 

In summary, to find out a consistent description for the ARuCVtype ruthenates, we study 
the temperature-dependent physical properties of recently synthesized cubic phase BaRu0 3 
by using the first principles LDA + DMFT approach. Judged from the quasiparticle band 
structure and integrated spectral functions of Ru-4<i states, the cubic phase BaRuOs is a 
weakly correlated Hund's metal. There exists local magnetic moment and the inverse local 
magnetic susceptibility obeys the Curie- Weiss law in the studied temperature regime. The 
low-frequency self-energy function, scattering rate, and optical conductivity of cubic phase 
BaRu0 3 show apparent NFL behaviors. It is argued that the Hund's rule coupling J plays 
an important role in the underlying physics in cubic phase BaRuOs and other perovskite 
ARUO3 compounds. 
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